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Abstract — A R&D project has been launched in 2009 to address 
fundamental methods in radiation transport simulation and 
revisit Geant4 kernel design to cope with new experimental 
requirements. The project focuses on simulation at different 
scales in the same experimental environment: this set of prob- 
lems requires new methods across the current boundaries of 
condensed-random-walk and discrete transport schemes. An ex- 
ploration is also foreseen about exploiting and extending already 
existing Geant4 features to apply Monte Carlo and deterministic 
transport methods in the same simulation environment. An 
overview of this new R&D associated with Geant4 is presented, 
together with the first developments in progress. 



I. Introduction 

GEANT4 JT], El is an object oriented toolkit for the sim- 
ulation of particle interactions with matter. It provides 
advanced functionality for all the domains typical of detector 
simulation: geometry and material modelling, description of 
particle properties, physics processes, tracking, event and run 
management, user interface and visualisation. 

Geant4 is nowadays a mature Monte Carlo system; its multi- 
disciplinary nature and its wide usage are demonstrated by the 
fact that its reference article [ 1 ] is the most cited publication 
0, El in the Nuclear Science and Technology category. 

Geant4 is the result of a R&D project (CERN RD44) 
carried out between 1994 and 1998. RD44 was launched at a 
time when the LEP experiments were running GEANT 3 as 
a well-established system, that had been refined throughout a 
decade of production service. RD44 investigated the adoption 
of the object oriented technology and C++ for a simulation 
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system replacing GEANT 3.21 [6|, and developed the first 
functional version of the Geant4 released at the end of 1998. 

The Geant4 releases following the first Geant4 one at the 
end of 1998 have added new functionality to the toolkit; never- 
theless, its architectural design and the fundamental concepts 
of Geant4 application domain have remained substantially 
unchanged since their original conception. 

New experimental requirements have emerged in the recent 
years, which challenge the conventional scope of major Monte 
Carlo transport codes like Geant4. Research in nanodosimetry, 
nanotechnology-based detectors, radiation effects on compo- 
nents in space and at high luminosity colliders, nuclear power, 
plasma physics etc. have shown the need of new methodolog- 
ical approaches to radiation transport simulation along with 
new physics functionality in Geant4. A common requirement 
has emerged in all such research domains, i.e. the ability to 
change the scale at which the problem is treated in a complex 
simulation environment. This requirement goes beyond the 
traditional issues of variance reduction, for which current 
Monte Carlo codes provide abundant tools and techniques. 

Significant technological developments both in software and 
computing hardware have also occurred since the RD44 phase. 
New software techniques are available nowadays, that were not 
yet established at the time when Geant4 was designed. 

II. Recent experimental evolutions and 

REQUIREMENTS 

Various experimental domains would profit of the capability 
of dealing with co-working transport schemes. 

Radiation effects at the nano-scale are important for the 
protection of electronic devices operating in various radiation 
environments. In realistic use cases small-scale systems are 
often embedded in larger scale ones: for instance, a component 
may operate within a collider experiment or on a satellite in 
space, cellular and sub-cellular aggregates in real biological 
systems exist in complex body structures etc. A simulation 
system capable of addressing different scales would be relevant 
to studies of the effects on components exposed to the fierce 
experimental environment of LHC (and future super-LHC). 

The capability of condensed and discrete Monte Carlo sim- 
ulation schemes in the same software environment is also crit- 
ical to experimental configurations involving nanotechnology- 
based detectors. While R&D in nanotechnology is actively 
pursued also in view of application to HEP detectors, it is not 
yet possible to simulate such detectors as standalone systems 
with Geant4, nor to evaluate their performance once they are 
embedded in a full-scale experimental set-up. 



Plasma physics requires addressing the concept of object 
state and behaviour mutation in relation to the environment: 
in this use case the mutability concerns both the physics pro- 
cesses and the particles involved. Astrophysics and studies for 
fusion-based nuclear power are just two relevant applications, 
which would profit of Geant4 applicability to this physics 
domain. 

The issue of co-existing condensed-random-walk and dis- 
crete schemes arises in another context of the simulation 
domain. It concerns the conceptually correct treatment of the 
atomic relaxation following the impact ionization produced 
by charged particles: since the cross section for produc- 
ing secondary electrons from ionization is subject to in- 
frared divergence, in the conventional condensed-random-walk 
schemes the interaction is treated in two different regimes 
of continuous energy loss along the step and of discrete 
X-ray production, with the consequent adoption of cutoffs. 
This scheme introduces an artificial dependency on cuts in 
the generation of PIXE (Particle Induced X-ray Emission), 
while atomic relaxation is intrinsically a discrete process; 
moreover, the current Geant4 scheme neglects the correlation 
between the X-ray spectrum of primary ionization and PIXE. 
A conceptual revision of the continuous energy loss and 
discrete scheme would be desirable in this physics domain too. 
Use cases affected by the current conceptual limitation to treat 
PIXE correctly in Geant4 involve multiple, multi-disciplinary 
domains: applications for material analysis from planetary sci- 
ence to cultural heritage, precise dosimetry, critical shielding 
optimization of X-ray telescopes etc. 

III. Main areas of research 

A R&D project, named NAN05, has been recently launched 
(7) to address fundamental methods in radiation transport 
simulation; it explores possible solutions to cope with the new 
experimental requirements mentioned in the previous section 
and evaluates whether and how they can be supported by 
Geant4 kernel design. It was initiated by a team at the Italian 
Institute of Nuclear Research (INFN) and currently gathers an 
international team of physicists and engineers with background 
in various disciplines. 

The main focus of the project lies in the simulation at differ- 
ent scales in the same experimental environment: this objective 
is associated with the research of transport methods across 
the current boundaries of condensed-random-walk and discrete 
transport schemes. An exploration is also foreseen about 
exploiting and extending already existing Geant4 features to 
apply Monte Carlo and deterministic transport methods in the 
same simulation environment. The Geant4 toolkit is the ideal 
playground for this research, thanks to the object oriented 
technology it adopted in the RD44 phase. 

Other issues have been identified along with the experience 
of Geant4 development and usage over the past 10 years, 
which would profit of the exploratory research in the kernel 
design motivated by the main objectives of the project: 

• Customization of physics modeling in a simulation ap- 
plication 

• Scattered and tangled concerns across the code 



> Facilities for physics verification and validation 

> Performance 

These topics are considered as supporting developments, 
which are instrumental to achieve the main goals of the project. 

The project adopts a software process model based on the 
Unified Process [8| framework. The software developments 
are motivated by concrete experimental applications, and sig- 
nificant effort is invested in the software design: these features 
of the project are well served by the Unified Process, which is 
use case driven and architecture-centric. The adopted software 
process framework involves an iterative and incremental life- 
cycle. 

IV. Co-working condensed and discrete 

SIMULATION METHODS 

Methods to model hard interactions of particles with matter 
constituents by means of an appropriate binary theory are well 
established: in this approach collisions are treated as binary 
processes, that is, either the target electrons are treated as free 
and at rest, or the influence of binding is accounted only in 
an approximated way. 

General-purpose Monte Carlo codes, like EGS ||9), iflOl . 
flU, FLUKA d, d, Geant4 and MCNP JH, fl3], flU, 
operate in this context. Their calculations of energy deposit 
distributions are based on condensed-random-walk schemes 
of particle transport. Charged particle tracks are divided into 
many steps, such that several interactions occur in a step; 
one energy loss and one deflection are calculated for each 
step. A further simplification consists in the adoption of the 
Continuous Slowing Down Approximation (CSDA), where 
the energy loss rate is determined by the stopping power. 
This approach is adequate as long as the discrete energy loss 
events treated are of magnitudes larger than electronic binding 
energies. 

Various specialized Monte Carlo codes, usually known as 
"track structure codes", have been developed for micro and 
nano-dosimetry calculations. They handle particle interactions 
with matter as discrete processes: all collisions are explicitly 
simulated as single-scattering interactions. This approach is 
suitable to studies where the precise structure of the energy 
deposit and of the secondary particle production associated 
with a track is essential; nevertheless, the detailed treatment 
of collisions down to very low energy results in a high 
computational demand. 

So far, simulation based on condensed-random-walk 
schemes and so-called track structure generation have been 
treated as distinct computational domains; this separation is 
due to the conceptual and technical difficulty of handling the 
two schemes in the same simulation environment. Achieving a 
conceptual approach and an architectural design where the two 
schemes can co-work would represent a significant progress 
in Monte Carlo simulation. 

Recently, a set of specialized processes for track structure 
simulation in liquid water has been designed and implemented 
in Geant4 [ 19 1; like their equivalents in dedicated Monte Carlo 
codes, they operate in the regime of discrete interactions. 
While the toolkit nature of Geant4 allows the co-existence of 



tools for simulation at different scales, the capability of these 
two schemes to effectively co-work in a multi-scale problem 
is still far from being established. 

V. Co-working Monte Carlo and deterministic 

SIMULATION METHODS 

Deterministic transport methods are widely used in various 
domains: radiotherapy treatment planning and calculations 
for nuclear reactors are just two examples. Their usage is 
motivated by the requirement of a fast simulation response 
in complex situations. 

Both Monte Carlo and deterministic transport methods have 
their own strong points and limitations. Research in this field 
to perform both transport models within the same simulation 
environment would be highly valuable: the capability of dif- 
ferent transport methods in the same simulation environment 
would simplify the geometrical and material modelling of the 
system under study, and would facilitate the analysis of the 
behaviour of the system itself. 

VI. Supporting research topics 

The complexity of the problem domain to be addressed 
requires the investigation of software techniques, capable to 
effectively support the conceptual objectives to be pursued. 

A. Generic programming techniques in physics simulation 
design 

Metaprogramming has emerged in the last few years as a 
powerful software technique. In C++ the template mechanism 
provides naturally a rich facility for metaprogramming; li- 
braries like Boost ifTTl and Loki lfT8l are nowadays available to 
support generic programming development. Metaprogramming 
presents several interesting advantages, which propose it as a 
worthy candidate for physics simulation design. 

This technique has not been exploited in Geant4 core yet: 
the evolution towards the C++ standard still in progress and 
the limited support available in C++ compilers at the time of 
the RD44 phase discouraged the consideration of templates 
for extensive use in Geant4 architectural design at that stage. 

A preliminary investigation of its applicability in a multi- 
platform simulation context has been carried out by one of the 
authors of this paper through the application of a policy-based 
class design [19| limited to a restricted physics sub-domain. 

B. Design for scattered and tangled concerns 

The problem domain of radiation transport simulation in- 
volves a number of concerns, which are common to multiple 
parts of the system, but whose code gets scattered across 
different parts; moreover, multiple concerns may be tangled in 
the same code. The capability of addressing scattered concerns 
by an effective design would results in leaner, more easily 
maintainable code: an optimization in this domain would be 
meaningful in a large-scale software system like Geant4. 

The object oriented technology lacks proper instruments 
to address the issue of scattering and tangling of concerns. 



Aspect oriented programming provides support for cross- 
cutting concerns (i.e. aspects) and for automatically propagat- 
ing appropriate points of execution in the code; nevertheless, 
this technology is not widely established yet, and language 
support is still relatively limited in C++. 

Two topics associated to concerns in Geant4 code are 
relevant to the research areas considered in this project: the 
issue of endowing objects - in particular, physics objects - 
of intrinsic verification and validation capabilities (more in 
general, of analysis capabilities), and dealing with secondary 
effects following a primary interaction (e.g. the relaxation of 
an excited atom). An exploratory study in these areas would 
be useful to evaluate the possibility of addressing concerns 
effectively in a simulation environment like Geant4. 

VII. Ongoing activities 

The activity currently in progress has an exploratory char- 
acter: it evaluates various problem domains to identify the 
issues to be addressed, the requirements in the associated 
experimental domains and candidate technologies. 

Regarding the main topic of research, the problem domain 
analysis has identified the concept of "mutability" as a main 
issue in the context of transition between co-working con- 
densed and discrete transport schemes. The current research 
in software design explores the introduction of the concept 
of "mutants" in the software design, and of "stimuli" capable 
of triggering mutations. Related concepts, like reversible and 
spontaneous mutation, are subject to investigation too. 

The introduction of the concept of mutability in physics- 
related objects requires the identification of their stable and 
mutable states and behaviour, and their fine-grained decompo- 
sition into parts capable of evolving, or remaining unchanged. 
Two pilot projects are in progress to explore the capability of 
policy-based class design to support this requirement in two 
large scale physics simulation domains: a general-purpose one 
ll20ll and a "track structure" one. 

Issues related to PIXE simulation are explored in another 
dedicated pilot project [21 J, [22|. 

In parallel, a project focussed on the simulation of radioac- 
tive decay [23 1 explores issues related to the collaboration be- 
tween electromagnetic and hadronic components of the design 
associated with the prototypes explored in the electromagnetic 
domain. 

VIII. Conclusion and outlook 

A R&D project is in progress to address the capability 
of dealing with multi-scale use cases in the same simulation 
environment based on Geant4: this requirement involves the 
capability of handling physics processes according to different 
transport schemes. The first developments associated with this 
project are described in greater detail in other contributions 
to the 2009 IEEE Nuclear Science Symposium Conference 
Record. M. Kuster and S. Hauf acknowledege support by the 
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Deutsches Zentrum fur Luft- und Raumfahrt - DLR under the 
grant number 50QR0902. 



Acknowledgment 

The authors thank Sergio Bertolucci (CERN), Simone Gi- 
ani (CERN), Vladimir Grichine (Lebedev Institute), Bernd 
Grosswendt (formerly PTB, retired), Alessandro Montanari 
(INFN Bologna), Andreas Pfeiffer (CERN), Reinhard Schulte 
(Loma Linda University), Manju Sudhakar (ISRO and INFN 
Genova) and Andrew Wroe (Loma Linda University) for 
helpful discussions. 

References 

[1] S. Agostinelli et al., "Geant4 - a simulation toolkit" Nucl. Instrum. Meth. 
A, vol. 506, no. 3, pp. 250-303, 2003. 

[2] J. Allison, et al., "Geant4 Development and Applications", IEEE Trans. 
Nucl. Sci., vol. 53, no. 1, pp. 270-278, 2006. 

[3] T. Basaglia, Z.W. Bell, P. Dressendorfer, A. Larkin, M.G. Pia, "Writing 
Software or Writing Scientific Articles?" IEEE Trans. Nucl. Sci., vol. 
55, no. 2, pp. 671-678, Apr. 2008. 

[4] M. G. Pia, T. Basaglia, Z.W. Bell, P. Dressendorfer, "Geant4 in scientific 
literature", in Conf. Rec. 2009 IEEE Nucl. Sci. Symp. 

[5] S. Giani, "GEANT4 : an object-oriented toolkit for simulation in HEP", 
CERN Report CERN-LHCC-98-044, Geneva, 1998. 

[6] "GEANT Detector Description and Simulation Tool", CERN Program 
Library Long Writeup W5013, 1995. 

[7] M. G. Pia et al., "R&D for co-working condensed and discrete transport 
methods in Geant4 kernel", in Proc. Int. Conf. on Mathematics, Compu- 
tational Methods and Reactor Physics (M&C 2009), on CD-ROM, Am. 
Nucl. Soc, La Grange Park, IL, 2009. 

[8] I. Jakobson, G . Booch, and J. Rumbaugh, "The Unified Development 
Process", Ed. Addison- Wesley, Boston, 1999. 

[9] W.R. Nelson, H. Hirayama, and D.W.O. Rogers, "The EGS4 code sys- 
tem", Report SLAC-265, Stanford Linear Accelerator Center, Stanford, 
CA, USA, 1985. 

[10] I. Kawrakow and D.W.O. Rogers, "The ENSnrc Code System: Monte 

Carlo Simulation of Electron and Photon Transport", NRCC Report 

PIRS-701, Sep. 2006. 
[11] H. Hirayama, et al., "The EGS5 code system", Report SLAC-R-730, 

Stanford Linear Accelerator Center, Stanford, CA, USA, 2006. 
[12] A. Fasso et al., "The physics models of FLUKA: status and recent 

developments", in Proc. Computing in High Energy and Nuclear Physics 

2003 Conference (CHEP2003), La Jolla, CA, USA, paper MOMT05. 
[13] A. Ferrari et al., "Fluka: a multi-particle transport code", Report CERN- 

2005-010, INFN/TC-05/11, SLAC-R-773, Geneva, Oct. 2005. 
[14] X-5 Monte Carlo Team, "MCNP - A General Monte Carlo N-Particle 

Transport Code, Version 5", Los Alamos National Laboratory Report 

LA-UR-03-1987, Apr. 2003, revised Mar. 2005. 
[15] RA. Forster et al., "MCNP Version 5", Nucl. Instrum. Meth. B, vol. 

213, pp. 82-86, 2004 
[16] J.S. Hendricks et al., "MCNPX, Version 26c", Los Alamos National 

Laboratory Report LA-UR-06-7991, Dec. 2006. 
[17] B. Karlsson, "Beyond the C++ Standard Library: An Introduction to 

Boost" YAddison-Wesley, 2005. 
[18] A. Alexandrescu, "Modern C++ Design", Ed.: Addison- Wesley, 2001. 
[19] S. Chauvie et al., "Geant4 physics processes for microdosimetry simu- 
lation: design foundation and implementation of the first set of models", 

IEEE Trans. Nucl. Set, vol. 54, no. 6, Dec. 2007. 
[20] M. Augelli et al., "Research in Geant4 electromagnetic physics design, 

and its effects on computational performance and quality assurance", in 

Conf. Rec. 2009 IEEE Nucl. Sci. Symp. 
[21] M. G. Pia et al., "PIXE simulation with Geant4", IEEE Trans. Nucl. 

Set, vol. 56, no. 6, Dec. 2009. 
[22] M. G. Pia et al., "Recent Developments on PIXE Simulation with 

Geant4" in Conf. Rec. 2009 IEEE Nucl. Sci. Symp. 
[23] S. Hauf et al., "Progress with Geant4-based Radioactive Decay Simula- 
tion Using the Example of Simbol-X and the IXO WFI and HTRS", in 

Conf. Rec. 2009 IEEE Nucl. Sci. Symp. 



